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We report on a study of topological properties of Fibonacci quasicrystals. Chern numbers which 
label the dense set of spectral gaps, are shown to be related to the underlying palindromic symmetry. 
Topological and spectral features are related to the two independent phases of the scattering matrix: 
the total phase shift describing the frequency spectrum and the chiral phase sensitive to topological 
features. Conveniently designed gap modes with spectral properties directly related to the Chern 
numbers allow to scan these phases. An effective topological Fabry-Perot cavity is presented. 
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Topological features of quasi-periodic (QP) structures 
have been studied under various guises both in physics 
and mathematics mil]. For instance, Chern integers la¬ 
bel the dense set of gaps in the energy (frequency) spec¬ 
trum. More generally, Chern numbers are known to play 
a role in problems where the underlying topology has 
been identified. Examples are provided by the quantum 
Hall effect, anomalies in Dirac systems and topological 
insulators. For each case, the origin of topological fea¬ 
tures is identified e.g., magnetic field, chirality or non¬ 
trivial band structures. However, in QP structures like 
the Fibonacci chain, the physical origin underlying the 
existence of Chern numbers has not yet been identified. 
The purpose of this letter is to address this point and to 
propose new expressions and possible measurements of 
Chern numbers using scattering data. We shall discuss 
the specific case of Fibonacci chains, but we expect our 
results to apply to a larger class of similar QP one-(or 
higher) dimensional structures. 

QP systems share similar spectral properties [IHZl,e.g. 
their singular continuous frequency spectrum has a frac¬ 
tal structure (Cantor set) [8], with a dense distribution 
of gaps characterised by the gap labelling theorem [T] . It 
provides a precise expression of the (normalised to unity) 
integrated density of states (IDOS) JV (sgap) at energies 
Sgap inside the gaps. For the Fibonacci chain, gaps are 
labeled by means of two integers (p, q) such that, 

^f{egap)=P + qT~^, ( 1 ) 

where r = (1 + is the golden mean. The integers 

q are Chern numbers [2] and p{q) keeps M {sgap) within 
[0,1]. These features which have been widely studied 
HHsuin] and recently measured using cavity polaritons 
HQ], are shown in Figj^a. 

Here, we consider finite Fibonacci chains and study 
the behaviour of conveniently emulated edge states (gap 
modes) driven by an angular variable (p (defined here¬ 
after) accounting for the important but not yet stud¬ 
ied palindromic or mirror symmetry of quasicrystalline 
chains. As we shall see, p whose role has been stud¬ 


ied in the Harper model or quantum Hall related setups 
[asimiii], allows to monitor this structural symmetry. 

To better understand the meaning of 0, it is instructive 
to review basic rules used to generate finite QP struc¬ 
tures. We study, for simplicity, sequences built using 
a two letters alphabet These letters may de¬ 

scribe different physical setups such as electromagnetic 
waves in a dielectric with refractive index modulation 
or a Schrodinger equation with a modulated potential 
HQ). We first consider the substitution algorithm (equiv¬ 
alent to the concatenation rules) [131 El- ^ substitu¬ 
tion a acts onto {A, B} according to cf{A) = AB and 
(j{B) = A. Successive applications of a generate a se¬ 
quence Sj = (7^{B) whose length is the Fibonacci number 
Fj>i = Fj-iPFj- 2 ^ with Fo =Fi = 1. The ratio Fjj^ijFj 
tends to the golden mean r in the limit j ^ oo. The 
corresponding sequence Soo becomes rigorously QP and 
invariant, i.e., self-similar under this iteration transfor¬ 
mation. The gap labeling relation 0 is established using 
this algorithm [T] (see for an elementary derivation). 
A second set of rules includes two alternative approaches. 
The first (HIH], relies on a characteristic function Xn tak¬ 
ing two possible values ±1 respectively identified to the 
letters {B,A}. Among possible choices and using a dis¬ 
crete analogue of the continuous Aubry-Andre model [5] , 
we consider the form. 


Xm = Sign 


cos ^27r mr ^-h — cos (tt r , (2) 


proposed in HD- The angular degree of freedom cp can 
be safely ignored for Soo but not, as we shall see, for 
the finite segment ^Ar = [xiX 2 * * 'Xa^] of The chain 
Sj generated using substitution is reproduced by setting 
P=Pf = 7tt~^, and N=Fj. We thus redefine ([^ using a 
shifted phase p = p — pF‘ 

Finally, the Cut&Project method (hereafter C&P) is 
used to investigate quasicrystals, tilings and dynamical 
systems [T3l[T4] (see Figj^. A QP chain is obtained from 
the Z^-lattice cut by a line A defined by y = T~^x-\-const. 
We denote the direction perpendicular to A and de- 
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FIG. 1: (color online) (a) Normalised IDOS for the Fibonacci 
chain [29] . The location of 8 selected gaps using their topolog¬ 
ical numbers [p, q] are indicated in accordance to 0 - (b) The 
corresponding n gap modes behaviour as a function of 

0 compared to (top), (c)-(d) p{k), obtained from 0 for 
'PN (c), and 'PnF N (d). 


fine an acceptance window by the band of width cen¬ 
tered at A. This realises the ”cut”. A C&P set is ob¬ 
tained by projecting the points inside O on A and 
along A^. The two possible distances along A between 
neighbouring projections, are denoted We thus 

generate ^Soo for which the choice of origin on A is irrel¬ 
evant. For finite chains it is not so, since the origin fixes 
the first letter and the iteration of the sequence. We 
note that C&P and characteristic function methods are 
related through the constant term in the equation for A, 
namely, y = T~^x—^. This allows for a useful interpreta¬ 
tion of 0 which accounts for the rearrangement of letters 
in a Fibonacci chain resulting from its translation along 
the ^-axis (FigJ^b). The 27r-periodic phase 0 is a struc¬ 
tural degree of freedom (a.k.a a phason). For a chain of 
length A, each value of 0 generates a different segment 
of 5'oo, and corresponds to a translation 6n = (r/27r)(/) 
along A. Monitoring (p induces a series of N identi¬ 
cal local structural changes equivalent to the inversion 
of a single 6-letters segment BAABAB ^ BABAAB 
(Figs0 a-b). These changes occur one at a time and are 
distributed according to a geometrical pattern displayed 
in Figj^c. This discreteness underlies the characteristic 
stepwise shapes visible in Figs{2b,[^d and[^ 

The combined description using p and its associated 
C&P interpretation allows to unveil an essential struc¬ 
tural symmetry of QP chains, namely their palindromic 
character (mirror symmetry). Sweeping p over a pe¬ 
riod, we observe that any chain of length N becomes 
perfectly symmetric for two specific values, (ppaiiN) = 
— {N + 2)7rr“^+m7r, where mGZ. Deviations from this 
palindromic symmetry is described using a structural pa- 





FIG. 2: (color online) (a)-(b) The G&P method for pp and 
0F —O.SStt. An ellipse indicates the origin of the chain. The 
resulting A = 10 sequences are in the insets. Shifting the 
band D along the vertical axis induces 2 structural changes 
corresponding to points entering and leaving D (red/green 
arrows), (c)-(d) Structure properties of m as a function 
of p. (c) Structural color plot. Black(white) regions describe 
B{A) letters. Identical structural changes (red circles, inset) 
occur when scanning p (e.g. green lines), (d) Behaviour of 
rf{p) (SiQ is indicated). The saturation at ? 7 ~ 0.76 is specific 
to Fibonacci through the [AA] pair occurrence, 2t —3 [15]. 

rameter defined in [0,1] by 

IXi+i(</’) - Xiv-i (</')! • (3) 

j=0 

The palindromic cycle for r]{p) is displayed in Figj^. 
This periodic occurrence of palindromic chains, as we 
shall see, is underlying the existence of Chern num¬ 
bers in the gap labeling relation 0. It is characteristic 
of C&P QP systems and of the corresponding Aubry- 
Andre-Harper tight binding model HE]. The structural 
phase p drives a ID structural symmetry cycle, and thus 
plays a role analogous to the magnetic field in the 2D 
quantum Hall effect. It is therefore natural to set the 
origin of p at the palindromic symmetry, i.e., to redefine 
(t> as (l)-(t>pai{N) in The observation of a ^-driven 
palindromic cycle is one of the main result of this letter. 

We now introduce the scattering formalism offering an 
elegant framework to study the spectral consequences of 
this structural symmetry m- We stress that although 
self-similarity and gap labeling 0 are spectral proper¬ 
ties of the infinite Soo (where p is irrelevant), they are 
equally relevant to finite chains (Figj^c). We recall that 
any quantum or wave system with a potential defined 
w.r.t a free part can be probed using scattering waves of 
wave vector k. Here, we consider a Fibonacci dielectric 
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FIG. 3: (color online) (a) The scattering setup. A finite chain 
is inserted in between perfect dielectric media supporting in¬ 
coming (zlA-r) and outgoing waves {ol^or). Displayed case: 
P = ^ N (1^ indicates the interface) (b) The chain ^n 
with a reflecting boundary condition (golden bar), is equiv¬ 
alent to the unfolded open structure P in (a), (c) Complex 
amplitudes corresponding to a wave incident respectively on 
the right or left side of the chain. 


medium with scattered electromagnetic waves [15] . With 
clear enough notations (FigJ^a), the (unitary) scattering 
matrix of the sequence F at is 

It is diagonalisable under the form diag{e^^'^ ^ We 

define the total phase shift 5{k) = (0i(/c) + (/)2(/c))/2. The 
transmission and reflexion amplitudes, t= = 

and V = |r|e*^T corresponding to left or right 
incoming waves are represented on Figj^c. The relation, 
det (k) = = —tjP^ implies that 6{k) = 0t{kyr7r/2. 

A simple relation exists between 6{k) and the density of 
modes p{k) [19]. For it reads 


showing that while the difference of the two reflected 
phases is the ^-dependent chiral phase a, their sum is 
the ^-independent total phase shift 6. 

The phase (a(0) is most conveniently investigated us¬ 
ing states other than ^-independent scattering states. We 
thus consider distinct choices of boundary conditions and 
study the resulting edge states whose spectral locations 
depend on these choices, e.g. closed (mirror) boundary 
condition, interface or local defect [21]. Since each choice 
leads to edge states of the same topological content [22] , 
without loss of generality, we consider a finite Fibonacci 
chain (of length N = Fj^ to properly fix the origin of (p) 
with reflecting boundary condition at one end (FigJ^b). 
To adapt this setup to the scattering approach, we unfold 
it as displayed in Fig^a, so that the resulting artificial 
palindrome P = ^n F n is an open scattering system of 
length 2N. Since P is not a bona fide Fibonacci chain, ad¬ 
ditional interface modes appear in the gaps of ^at at val¬ 
ues km{<p) whose behaviour depends solely on the Chern 
number q of each gap as assigned by relation 0- Modes 
cross the gaps in a direction set by sign{q) and a crossing 
period of 7r/\q\, as shown in Figsj^b,d [23]. When ^at 
is palindromic (^at = ^Ar), P reduces to F at with peri¬ 
odic boundary conditions instead of a mirror. Therefore, 
no gap modes exist (Figj^b,d). 

These observations become more quantitative by ex¬ 
pressing the scattering data of the open structure P in 
terms of those of ^at- The respective scattering matrices 
^ of the chains P, ^at and ^at, are related by, 

detPp = (detP^) (detP^) , (7) 


l>{k) - Po{k) = ^ I- A l„detS,(t) = . (5) 

where po{k) is the density of modes of the free system, 
i.e. withuot ^ at- This relation is used to obtain the 
spectrum of Figj^c. Since p{k) is independent of the in¬ 
cident wave direction, the phases S and dt do not dis¬ 
criminate between the sequence ^at and its reversed 
F N = [Xat Xat-i* •'Xi], and therefore they are not sen¬ 
sitive to the symmetry of the sequence (palindromic or 
not). This insensitivity, in addition to the fact that shifts 
in (p are only shifts along Poo, clarifies why these phases, 
describing bulk properties, are independent of p [20] . 
This is not so for the reflexion amplitudes, related by 
V = ^ where a= 0^—0^ is a chiral phase which, 

by definition, vanishes identically for a palindromic ^ at 
and is finite otherwise. Within the gaps \a{p,k) \ is a 
function of essentially insensitive to k (FigsJ^,e), thus 
making it a spectral counterpart of p{p). Finally, the uni- 
tarity condition, tr^*+Vt* = 0, of in Q, provides the 
relation, 

20t + 7r = 0 ^ F 0 ^ = 25 ^ (6) 


where = (1 —^*)/( 2 : —1) is a Fabry-Perot factor with 

z = = |rp ^ F) = |rp [15]. Since 

—y i — 

by construction, 0^ = the cavity phase Ocav^ a 

denomination we shall justify later on, becomes 

dcav = ^ F 0^ = 20^. ( 8 ) 

Moreover, inserting the relation det = det 
into 0> leads to 5p = 25 F (pf 2^ 5p being the phase 
shift of Pp. A condition for the appearance of a gap 
mode km in the spectrum of P is 5p{km) — ‘^5{km)=7TTn^ 
where m G Z or equivalently p = 27rm [24] . In the limit 
|rp ^ 1 [15], we obtain for the Fabry- 

Perot factor, so that the resonance condition on p be¬ 
comes Ocav {km) = ‘2.Tim. Since, according to (|^, Ocav is a 
phase of P^, the gap modes km{P) of ^at bounded by 
a mirror at one end, can be obtained from the scatter¬ 
ing data of the corresponding open chain [23]. Knowing 
dcav (0, k) allows to determine precisely the spectral lo¬ 
cation km{P) of the gap modes during a palindromic cycle 
and to characterise their crossing direction and period¬ 
icity as a function of the corresponding Chern number q 
(Figs|^c,d,g,h). This quantitative description is the basis 
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FIG. 4: (color online) The phases aq and Ocav for two gaps of m, with —1,2. (a) Color plot of |(a_i(0,/c)| within the 
gap. (b) Comparison between the midgap value, |(a-i(0)|, and ?7(0). |cKq(0)| closely follows the palindromic cycle, and is a 
7r/|g|-periodic spectral counterpart of the 7r-periodic ?7(0) [30] (c) Color plot of 0cav{(t>^k) within the gap q — —l. Light green 
areas correspond to the resonance condition 0cav = ‘^7rm. (d) Crossover of the gap mode km {(!>)■ Relative spectral location (red) 
is compared to the resonance condition (blue) and to ?7(0) (green), (e)-(h) Similar plots for q = 2. 


of the qualitative analysis displayed in Figsj2b,d. Using 
and the definition of the chiral phase a we obtain, 

0cav{4>,Q,k) =2S{k) +aq{(p) , (9) 


which disentangles Ocav into (5, insensitive to and a k- 
independent aq [25| , which displays the topological prop¬ 
erties, namely the 7r/\q\ period in (j). These properties do 
not depend on the structural modulation strength [T5] . 

As a straightforward consequence of the definition of 
Ocavi its winding number 

W {Ocav) = ^ ^^ (10) 

defined in a piecewise manner, provides a direct deter¬ 
mination of the corresponding Chern number q [26]. In 
addition, the ^-independence of J, together with (|^, lead 
to W {aq) = W {Ocav) = ‘^Q for the chiral phase wind¬ 
ing number. These expressions are especially interesting 
since they relate the Chern numbers labeling the spectral 
gaps of the infinite chain Soo to the scattering matrix of 
finite chains. It is thus possible to directly deduce Chern 
numbers from a scattering experiment. 

Finally, we provide another interpretation of the res¬ 
onance condition Ocav {km) = 27rm. The interface be¬ 
tween ^AT and ^TV defines a zero length effective cavity 
between two topological (g-dependent) and frequency- 
dependent highly reflective mirrors. The accumulated 
phase shift of this effective cavity is precisely the cavity 
phase 0cav=~^^ (Figj^c), thus justifying its name. 
We define a frequency-dependent effective cavity length. 


C{(l),q,k) 


m ^cav k) 
4 TT 


( 11 ) 


where \{k) = 271/k is the wavelength. Resonant modes 
occur at gap frequencies satisfying the usual Fabry-Perot 
condition 2>C(/c^)/A(/c^) = m G Z, i.e. the previous res¬ 
onance condition. This interpretation can be extended 


to situations where, in addition, a geometric cavity of 
length L is inserted between the two mirrors so that the 
effective length becomes C (0, g, km) = L^\{k) Ocav/ (^tt), 
thus leading to a new set of topological gap modes m- 

In summary, we have shown that topological proper¬ 
ties of Fibonacci QP chains encoded in Chern numbers 
labelling the dense set of spectral gaps, are related to the 
existence of an underlying palindromic symmetry driven 
by a structural gauge phase (j). Using a scattering ap¬ 
proach, we have shown that these Chern numbers can 
be obtained either by following the behaviour of con¬ 
veniently generated gap modes or by measuring scat¬ 
tering phases. Equivalently, finite length QP Fibonacci 
chains act as topological mirrors so that cavities defined 
by two such mirrors enclose modes whose properties are 
determined by topological features. This observation is 
of interest in different contexts such as Casimir physics. 
These results are readily applicable to the Aubry-Andre- 
Harper model. Finally, we expect some of the consider¬ 
ations presented in this letter to apply to QP systems 
in higher dimensions such as Penrose or Rauzy tilings 

[ZzlEH]. 
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SUPPLEMENTAL MATERIAL 

Notation and details on the scattering approach 

The scattering approach used in the letter is based 
on [17]. In this section we provide details, specific to the 
calculation of phases of the scattering matrix and spectra 
of the structures P, and F at- 

Each of these structures is considered as a ID dielectric 
with a uniform refractive index everywhere except for a 
finite spatial region, where this index is modulated. The 
two possible types of layers, {A, P}, are characterised 
by their refractive indices ub} and widths ds} 
such that riAdA = The latter condition 

is the quarter wavelength layer setup which is known to 
enhances the spectral effect of quasi-periodicity. 

For a transverse-electric electromagnetic wave prop- 
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FIG. 5: (color online) Setup and wave propagation notations 
for the dielectric scattering geometry problem, {il^ol} = 

I denote the electric field amplitudes of incom¬ 
ing and outgoing waves at the left boundary, respectively. 
Similarly, {zi?,Oi?} are defined outside the 

right structure boundary. 


agating in the modulation direction z, the scattered 
quantities are electric field amplitudes, E{z)^ such that 
E{z,t) = E{z)e~^^'^. The setup and the corresponding 
wave propagation notations are given in Figj^ 

We obtain the scattering matrix defined in the letter 
by (4), through the algebraically related transfer matrix, 
M(F), which is defined by 


Or 



1 (t/t* 

t 1 ) 



The transfer matrix of T at is calculated using its multi¬ 
plicative properties through 

As stated in the letter, the normalised IDOS of a struc¬ 
ture of length N is given hy N'{k) = This implies 

that for every new mode, the total phase shift S{k) expe¬ 
riences a discrete jump of size tt. This understanding will 
be used in subsection [] to derive the gap modes resonant 
condition. 


responding S^. To that purpose, we make use of the 
multiplicative properties of the transfer matrix. 


Mp 


1 /tp* ^pjtp 
-rp/tp 1/tp 


= . 


(14) 


Using (12), and recalling that = t, Mr becomes 


Mp 


1 / t/t* \ (t/t* 


(15) 


Using ( 6 ) of the letter, with 0 and here, the total 
phase shift (5p, may be extracted from^'^^ = —tp/tp = 
/^P• Noting that and , we 

arrive, after some straightforward algebra, to (7) of the 
letter. 


Derivation of the resonant eondition Ocav {km) — 27rm 

The appearance of new modes in the spectrum of the 
structure P is understood from the total phase shifts dif¬ 
ference Sp—26. This results from the fact that in the ideal 
case where |rp = 1 , i.e. when 6 = const^ the appearance 
of every new mode yields a discrete jump of size ir in 6p 
(see section I above). Therefore, when |rp = 1, the spec¬ 
tral locations of new gap modes may be expressed by the 
resonant condition 6p — 26 = 7 rm, which according to (7) 
of the letter yields {km) = 27rm. 

Next, we find the resonant condition for any |rp 7 ^ 0. 
Using (7) of the letter, one obtains 


|r|^ cos {Ocav + 2 ) = COS (|) , (16) 

or alternatively. 


scattering analysis for the artificial palindrome 

In this section we describe (A) the derivation of the 
phase ip defined in (7) of the letter, (B) the derivation 
of the resonant condition Ocav {km) = 27rm, and (C) the 
dependence of this condition upon dielectric contrast. 


Derivation of the phase p 

To derive the phase p appearing in (7) of the letter, 
one needs to relate the total phase shift of the scattering 
matrix of the structure P, 


to the total phase shift of the substructures scattering 
matrices defined in (4) of the letter, and the cor¬ 




From ( |16[ ) , 

if {0. 


one can see 


rr = l) = -Sc 


that for perfect reflectance 
which yields the condition 
for the appearance of new modes Ocav {km) = 27rm, as in 
the letter. Using ( pT[ ), FigJ^ shows the effect of varying 
|rp on the relation between p, and Ocav- ^ no longer cov¬ 
ers the interval [0, 27r], and consequently never takes the 
perfect reflection resonant values (while Ocav does). In¬ 
stead, it becomes increasingly smeared as |rp decreases. 
This shows that in general, Ocav is more suitable than p 
for identifying the emergence of new modes in structures 
with imperfect reflection. 

Figure shows the same conclusion for the dielec¬ 
tric structure [27]. The behaviour of Ocav = 2 6 > and 
p = 2Sp—4:S is examined for two different spectral gaps, 
characterised by different values of |rp. Indeed, the wind¬ 
ing range of Ocav is unaffected by |rp, while p ceases to 
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FIG. 6: (color online) (fi (9cav,\r\^) as calculated from 
The relation if = —Ocav holds only for perfect reflectance. 




0.365 0.37 0.375 0.38 0.385 0.39 0.395 
2nd/X 


FIG. 7: (color online) Behaviour of the phases Ocav {k) and 
f {Ocav {k ), \r{k)\‘^) as a function of k for the dielectric struc¬ 
ture [27], for different spectral gaps (gap edges are indicated 
by red bars), (a) Gap q — 2, with |r|^~0.25. (b) Gap q — —l, 
with |r|^?^0.98. 



FIG. 8: (color online) The spectral location of gap modes 
km{4>) as a function of contrast for the ^ — 1 gap of the 

structure P [27] at (a) 15% and (b) 80% dielectric contrast. 


in the letter), and of the Ocav and a are found to be un¬ 
changed by the dielectric contrast. However, the stepwise 
path taken by the gap modes while crossing the gaps, 
and the corresponding stepwise behaviour of the phases 
change significantly with the contrast. The discrete spec¬ 
tral changes in the phases and in the gap mode spec¬ 
tral location are not of equal magnitude as some changes 
are larger than others (this can be seen for instance in 
Figs.5.a,d of the letter). When the contrast is increased, 
the larger spectral jumps increase, while smaller spectral 
jumps decrease. For too large dielectric contrast, the en¬ 
tire crossover is accomplished by one spectral jump, thus 
unmeasurable. This trend can be seen in Figj^ 


Effective Fabry-Perot model with a cavity 


cover the interval [0,27r] (even for |rp = 0.98, Fig0b). 
Thus, for any value of |rp, the condition Ocav {km) = ^irm 
can be used to calculate the gap mode frequencies for P. 


Modulation strength dependence 

In dielectric structures, modulation strength is given 
by the dielectrc contrast, {nhigh — '^iow)/'^iow In struc¬ 
tures with spectral gaps, such as the Fibonacci chain, 
increasing the dielectric contrast results in an increase in 
gap widths and in the reflectance at gap frequencies. It is 
thus natural to ask how does the dielectric contrast affect 
the spectral topological properties discussed in the let¬ 
ter. The 7r/\q\ period of the gap mode crossover (Fig.lb 


In this section, we generalise the effective Fabry-Perot 
interpretation to include the case of a genuine cavity in¬ 
serted between the topological mirrors discussed in the 
letter. In the standard description, a Fabry-Perot cav¬ 
ity is defined by two mirrors separated by the length L,, 
and a discrete spectrum obtained by the resonant con¬ 
dition 2LlXm = Tn, may be viewed as the winding of a 
(frequency dependent) cavity phase. Op = dirL/X. This 
leads to a resonance condition Op {km) = 27rm. In the ef¬ 
fective Fabry-Perot model for the structure P, waves are 
also localized inside a finite region (the interface region) 
by means of distributed feedback, without a geometri¬ 
cal cavity, but only the sum of the two mirrors reflected 
phase shifts, Ocav Substituting Ophy Ocav yields the defi¬ 
nition of the virtual (frequency dependent) cavity length 
























































given by (11) of the letter. Figures]^ a andj^b display 
the predictions of the resonant condition for palindromic 
and non-palindromic values of respectively. 

As discussed in the letter, the effective Fabry-Perot 
model may be generalised to the case where in addition 
to the previous case, a genuine cavity of length L with a 
uniform refractive index, is inserted in between the topo¬ 
logical mirrors. This genuine cavity should be added to 
the effective one. In this case, the resonances may be 
found in the cavity phase picture as 


^total (^m) — ^cav (^m) T (^m) — 2i7r7Tl^ (1^) 

or alternatively in the cavity length picture given by 


‘^^total (^m) /^ (^m) ^5 

with Lfotai = L -\- C. As in the standard Fabry-Perot 
cavity case, the increase in the effective cavity length 
reduces the resonant mode separation. Since the gap 
widths did not change, multi-resonant-modes appear now 
in the gaps (FigJ^c). 


The gap labeling theorem 

In the following section, we present an elementary 
derivation of (1) in the letter, following the gap label¬ 
ing theorem [1] which applies to quasicrystals which can 
be constructed by substitution rules. Considering a 2 let¬ 
ter alphabet {A, 5}, and the substitutions {(7(A), 

(given in the letter), we first we define the one and two 
letters occurrence matrices, M and M 2 respectively. The 
matrix element contains the occurrence of the 
letter in the substitution. For the Fibonacci chain, it 
is 


(a) 


c/D 

O 2 
Q 


— In(DOS) 

o Fabry-Perot resonance 




0.35 0.4 0.45 0.5 0.55 0.6 0.65 

2nd/^ 


(b) 


C/D 

O 2 
Q 

^ 0 


— In(DOS) 

o Fabry-Perot resonance 




0.35 0.4 0.45 0.5 0.55 0.6 0.65 

2nd/;^ 


(b) 


O 2- 



^ —In(DOS) ^ 

0 Fabry-Perot resonance 

T' 

r( 

Ja ^^^l'. 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 

IndlX 


FIG. 9: (color online) (a)-(b) Density of states for the struc¬ 
ture P [27] (blue line). Red circles mark the Fabry-Perot 
resonant wave-vector values km, obeying Ocav (km) = 27r?n. 
(a) 0 = 0, yielding a palindromic structure, (b) A non- 
palindromic structure (0 = I.Ttt). (c) The genuine cavity 
setup: the substructures n, and ^n are separated by a 
uniform region with a refractive index n = ua, and width 
dfree^SOnd. The density of states (blue lines) and the reso¬ 
nant condition prediction (red circles) are given for 0 = I.Ttt. 


(72(AA) = (72(A5) = {AB){BA);a2{BA) = (AA), (21) 


(20) 

The matrix element contains the occurrence of the 
2-letter combination (doublet) in the 2-letter sub¬ 
stitution. For the Fibonacci chain, the 3 doublets alpha¬ 
bet is {AA, AB, BA}, and the corresponding substitu¬ 
tions are {cf{AA) = ABAB, a{AB) = ABA, a{BA) = 
AAB}. No BB doublets occur. The occurrences of the 
doublets in the 2-letter substitutions are obtained as fol¬ 
lows [1]. The numbers of doublets to be deduced (by 
shifting a 2-letter acceptance window from left to right) 
is equal to the number of letters in the substitution for the 
single first letter of the doublet. For a{AA) and a{AB) 
we should deduce 2 doublets (as the length of (j{A) = AB 
is 2), and for a{BA) we should deduce 1 doublet (as the 
length of (t{B) = A is 1). Finally, 


which leads to the occurrence matrix 


M 2 = 1 1 0 . (22) 

yi 1 B 

The next step is to diagonalize M and M 2 , and normal¬ 
ize the eigenvectors corresponding to the highest mutual 
eigenvalue, so that the sum of their elements is unity. The 
eigenvalues of M and M 2 are {1—T,r} and {1—T,r,0}, re¬ 
spectively. Thus, r is the mutual highest eigenvalue, and 
the corresponding normalized eigenvectors Vi, V 2 are 



( 23 ) 
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Using the identity r ^=t — 1, (12) becomes 


Now, Mgap^ the IDOS at possible gaps is equal to [1] 


/ 2r-3 
Uo = 2 - r 


The elements {^2 ) represent the occurrence of 
the letter (doublet). From Vi we see that the let¬ 
ter A occurs more than the letter B by 61.8% : 38.2%. 
From V 2 we see that the AB^ and BA doublets have 
an equal occurrence of ^ 38.2% each, and the A A dou¬ 
blet has ^ 23.6% occurrence. This clarifies the asymp¬ 
totic value of the plateau in 77 ( 0 ) (Fig.2.d in the let¬ 
ter). For the ^23.6% ^4^4 doublet occurrence in a given 
structure length there can never be a corresponding BB 
doublet preventing 7^(0) from approaching unity (anti- 
palindrome) . 

The normalized IDOS in the gaps can be calculated as 
follows [1]. Let Z[x] be the set of finite polynomials of 
degree N with integer coefficients. 






Ngap = { aP p-1) I P p-1) e Z } mod. 1, (26) 

where a G {Ui, U2}, and ^ is the highest common eigen¬ 
value of the occurrence matrices (i.e. 'd = r for the Fi¬ 
bonacci chain). Using (12), the set appearing in (15) for 
the Fibonacci chain is reduced to 


M -la —-a fU 
^gap - pi ! 93 


where ni^ 2 , 3 £N> and Since r ^ = 1 —t the 

set can be further compacted into 

•Vgap = {qT~^,q e z} mod. l, ( 28 ) 

where q e Z. This form leads to (1) of the letter. 



